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摘要 : 高 密 区 核 物质 状态 方程 的 研究 是 当前 核 物 理 研究 的 前 沿 热 点 问题 之 一 。 本 文 利用 极端 相对 论 量子 分 子 
动力 学 (UrQMD) 模型 ， 以 两 x 介子 Hanbury-Brown-Twiss (HBT) 关联 为 例 讨论 了 2-5 倍 饱 和 密度 (oo) 区 核 物 
质 状 态 方程 对 两 粒子 关联 效应 的 影响 。 通 过 使 用 不 同 的 核 物质 状态 方程 ， 展 示 了 密度 依赖 的 势 相互 作用 以 及 
TT 相 变 对 两 zx 介子 HBT 关联 和 zx 介子 发 射 源 时 空 必 性 的 影响 。 结 果 显示 在 4 ~ 5po 以 下 ,zt 介子 发 射 源 的 HBT 
A 半径 及 参数 敏感 于 核 物 质 状 态 方程 的 软 硬 。 通 过 与 实验 数据 的 比较 ， 现 有 的 HBT 实验 数据 排除 了 ~ 4po 以 下 
发 生 一 阶 相 变 的 可 能 性 ， 并 支持 一 个 在 低 密 区 CS 4p0) 表现 偏 硬 ， 在 高 密 区 由 于 相 变 而 逐渐 软化 的 核 物 质 状 态 
方程 。 研究 结果 强调 了 zz 介子 发 射 源 的 HBT 半径 及 参数 敏感 于 核 物 质 状 态 方程 的 软 硬 ， 可 用 于 理解 和 约束 高 
重子 密度 区 的 核 物质 状态 方程 。 
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Q 1 5 引言 变 弹 靶 核 的 质量 、 碰撞 能 量 和 碰撞 中 心 度 , 人 们 可 以 通过 
SN HICs 产生 不 同 密度 、 温 度 等 条 件 下 的 核 物 质 。 为 了 研究 
a +A YZ aya ET IDR yu H E 
D queer Ro AME QCD 相 图 结构 和 确定 CEP 在 QCD 相 图 中 的 位 置 ， 在 过 
SC 3 了 (Quantm Cromo ynamcsy 去 的 20 多 年 间 , 世界 上 众多 重 离子 对 撞 机 上 进行 了 一 系 
4 QCD) 相 图 和 相 结 构 的 确定 已 成 为 中 高 能 核 物 理 、 核 天 列 的 实验 ee A i 
CU M an Scanian tees : 9 实验， 如 位 于 美国 布鲁克 海 文 国家 实验 室 的 相对 论 
EC ” 体 物理 等 领域 的 前 沿 热点 课题 中 ,从 第 一 性 原理 出 发 的 "m pm ; 
$= 和 M 重 离子 对 撞 机 (Relativistic Heavy Ion Collider, RHIC) 上 进 
e ” 格 点 QCD 计算 表明 ， 在 零 重子 化 学 势 Up =0MeV) 下 z= 束 流 能 量 扫描 实验 ( 二 3 ~ 200 GeV J Au--A 
‘OBE SSIES ATI (Quark-Gluon-Plasma, QGP) ee ee A | 
EGP SIL ei Te iie 碰撞 ， 其 中 7.7 GeV 以 下 为 国定 靶 模 式 )a0。 通 过 对 多 
相 之 间 的 相 变 是 连续 过 渡 , 随 着 重子 化 学 势 的 升 高 , 格 点 0 oe poe PEL SERITUR EZ Z HOSP 

se ui 、 a M 个 敏感 于 关联 长 度 的 观测 量 进 行 分 析 ， 人 们 将 感 兴趣 的 
展 计算 83。 基 于 QCD 的 有 效 模型 计算 表明 在 较 大 重 Ro ele se. 人 na aa 

范围 内 5122。 此 外 ， 世 界 上 一 大 批 拟 建 以 及 新 建 的 加 


ay Mg bs tne HH as Al etea Man 而 一 

anu a non PISS 速 器 装置 也 拟 开展 质心 系 能 量 为 几 个 GeV 能 区 的 重 离 
yf HJI AM IKA AITA (Critical End Point, a eros As DJ TE E Er n ATLAS E 
CEP, fA, RAKERA PARUM muA ORTER E AAEE TIF”. 
置 仍 具有 很 大 的 不 确定 性 ， 存 在 明显 的 模型 依赖 性 eal。 理论 上 ， 输 运 模型 作为 抽取 核 物质 属性 的 重要 工具 


重 离子 倍 撞 (Heavy Ion Collisions, HICs) 实验 为 人 们 之 一 ， 经 过 多 年 的 不 断 发 展 和 完善 ， 已 经 形成 了 多 种 基 
研究 极端 条 件 下 核 物 质 的 性 质 提 供 了 条 件 10。 通 过 改 ” 村 不 同 理论 假设 和 处 理 的 用 于 描述 从 低能 到 相对 论 能 区 
HICs 的 输 运 模 型 n06131, 为 了 理解 和 降低 理论 系统 误差 以 
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在 抽取 核 物 质 属性 的 过 程 中 ， 实 验 测量 和 理论 计算 2 UrQMD 模型 及 HBT 简介 
结果 的 对 比 需要 基于 同一 观测 量 ， 因 此 构建 和 选取 敏感 "nM 
于 核 物质 状态 方程 Equation of State, Eos) 的 观测 量 尤 为 ”2.1 UrQMD 模型 简介 
重要 。 人 们 构建 并 提出 了 多 种 可 观测 量 用 以 探索 QCD fH 极端 相对 论 量子 分 子 动力 学 〈Ultrarelativistic Quan- 
边界 和 寻找 CEP 的 位 置 ， 如 集体 流 、 守 恒 荷 的 高 阶 累积 tum Molecular Dynamics, UrQMD) 模型 起 源 于 上 世纪 九 
和 矩 、 轻 核 的 产 额 比 、 间 歇 分 析 、 强 子 的 Hanbury-Brown- 十 年 代 ， 经 过 多 个 研究 组 持续 不 断 的 更 新 ， 现 已 被 广 
Twiss (HBT) 关联 6 9。 本 文 将 重点 关注 两 zx 介子 的 HBT —— 泛 的 用 于 模拟 从 费 米 能 区 到 LHC 能 区 相对 论 重 离子 碰 
关联 ， 它 可 以 用 于 揭示 HICs 中 冻结 时 刻 x 介子 发 射 源 $8 006371. E UrQMD 模型 框架 下 ， 具 有 一 定 宽度 的 高 斯 
的 时 空 属 性 和 x 介子 之 间 的 动量 相关 性 19。 关于 强 子 波 包 所 表示 的 强 子 将 在 平均 场 和 碰撞 项 的 共同 作用 下 进 
HBT 关联 的 详细 历史 和 发 展 介 绍 ， 请 感 兴趣 的 读者 参阅 行 演化 ， 第 i 个 粒子 的 坐标 r, 和 动量 了 的 时 空 演 化 遵循 
文献 L931,。 研究 表明 粒子 发 射 源 的 半径 (也 被 称 为 HBT ”哈密 顿 运动 方程 : 


半径 ) 敏感 于 EoS， 在 临界 点 附近 时 由 于 相 变 的 发 生 状 上 òH a 
态 方程 将 会 变 得 非常 软 ， 导 致 系统 膨胀 放 缓 效应 的 出 现 ， a Pr 


而 这 一 效应 可 以 通过 x 介子 发 射 源 的 HBT 半径 之 比 或 Hp, = YH, 是 系统 的 总 哈密 顿 量 ， 第 i 个 粒子 的 
者 平方 差 的 激发 函数 的 非 单调 (出 现 极 大 值 ) 行为 得 以 ” 哈密 顿 量 由 动能 和 势能 构成 H, = EF" 4. U,, 
TS Pel, 默认 版 本 的 UrQMD 模型 中 势能 包含 两 体 Skyrme 
势能 、 三 体 Skyrme 势能 、Yukawa 势能 、Pauli 势能 和 
基于 RHIC 束 流 能 量 扫描 实验 数据 , STAR 实验 组 对 ”库仑 势能 B34。 在 几 个 GeV. 以 上 能 区 的 重 离子 碰撞 中 
两 x 介子 HBT 关联 进行 了 分 析 , 结果 表明 在 VSNN m 20 ”Yukawa 势能 、Pauli 势能 对 整个 动力 学 演化 过 程 影响 较 
GeV Ik}, 2 介子 发 射 源 的 HBT EEE AEP, ”小 B89, 因此 在 本 文 的 研究 中 被 忽略 。 密度 依赖 的 Skyrme 
然而 通过 对 比 现 有 的 HBT 实验 数据 可 以 发 现 x 介子 发 ”势能 Uskyme 可 以 通过 密度 依赖 的 单 粒 子 势 得 到 并 表达 
射 源 的 HBT 半径 之 比 在 VsNN ~ 4 GeV 时 也 存在 着 极为 : 


大 值 ， 但 具有 极 大 的 不 确定 度 必 -29]。 理 论 上 ， 目 前 基 了 ( pp ) n (2 |: Q) 
不 同 模型 的 强 子 HBT 关联 的 研究 主要 集中 在 高 能 重 离 um Po Po 
子 碰撞 领域 中 I。 此 外 ， 对 于 基于 同一 实验 数据 或 者 输 。 ”其 中 po = 0.16 fm-3 WEI, p, 为 重子 密度 。 
运 模型 得 出 的 不 同 观测 量 的 碰撞 能 量 依赖 关系 中 ， 其 非 ” 通过 选取 不 同 的 势能 参数 可 以 给 出 不 同 的 状态 方程 。 鉴 
单调 行为 出 现 的 能 量 区 间 也 存在 一 定 的 差异 。 因 此 ,发 。 于 目前 人 们 对 高 密 区 核 物 质 状态 方程 的 软 硬 尚未 取得 统 
= 展 一 套 能 够 自 洽 地 描述 大 入 射 能 量 范 围 内 的 HICs, 并 可 一 的 认 知 02， 本 文采 用 不 可 压缩 系数 Ky = 200 MeV 和 
^ ”以 描述 多 个 观测 量 并 直接 与 实验 数据 对 比 的 输 运 模 型 显 Ky = 380 MeV 两 种 较 软 和 较 硬 的 状态 方程 , 相应 的 参数 
” ”得 十 分 重要 。 最 新 升级 的 UrQMD 模型 可 以 引入 任意 的 。” 如 表 1 中 所 示 。 需要 指出 的 是 ， 本文 暂时 没有 在 平均 场 中 
密度 依赖 的 状态 方程 下 ， 这 为 自治 地 研究 密度 依赖 。 考虑 x 介子 -核子 相互 作用 势 , 因为 致密 核 物质 中 的 x 介 
的 势 相互 作用 、 不 同 密度 条 件 下 的 相 变 对 不 同 观 测量 的 。 子 .核子 相互 作用 势 目前 为 止 还 没有 得 到 很 好 的 理解 ， 尚 


影响 提供 了 可 能 。 文献 中 基于 UrQMD 模型 展示 了 碰撞 。 存在 一 定 的 争议 B9421, 而且, CERAM SAAT Pe A, 
能 量 为 VsNN = 2 ~ 10 GeV HUSA PRET AAS 。 通过 在 UrQMD 模型 中 引入 不 同 强度 的 xN 相互 作用 势 ， 
依赖 的 势 相互 作用 以 及 相 变 对 集体 流 的 影响 ， 其 结果 排 EM x 介子 的 集体 流 及 其 中 心 度 、 快 度 依赖 的 实验 数据 
除了 4 倍 饱 和 密度 po 以 下 发 生 相 变 的 可 能 性 。 因 此 , 我 。 支持 一 个 较 弱 的 xN 相互 作用 势 。 


们 期 待 密度 依赖 的 势 相互 作用 也 将 对 重 离子 碰撞 中 粒子 表 1 Skyrme 势能 参数 
关联 效应 带 来 明显 的 影响 ， 从 而 澄清 ~ GeV 能 区 z 介子 一 一 一 一 一 一 
发 射 源 HBT 参数 的 能 量 依赖 谜团 并 确定 该 能 区 核 物 质 参数 。 Hard EoS Soft EoS 
状态 方程 的 密度 依赖 行为 。 a -124 -356 
f 71 303 
y 2.0 1.17 


本 文 框架 如 下 : 在 第 二 节 中 对 UrQMD 模型 和 两 x 
介子 HBT 关联 进行 简要 描述 ， 在 第 三 节 中 展示 核 物 质 除了 上 述 基 于 Skyrme 模型 给 出 的 EoS 外 ， 基 于 考 
状态 方程 对 x 介子 发 射 源 的 HBT 半径 的 影响 ; 最 后 , (E — 虑 不 同 强 相 互 作用 相 变 模式 的 手 征 平均 场 (Chiral Mean 
第 四 节 中 给 出 总 结 与 展望 。 Filed, CMF) 模型 的 EoS 也 被 引入 到 了 UrQMD 模型 


第 00 HA 李鹏 程 等 : CNPC2023: 利用 两 x 介子 HBT 关联 约束 重 离子 碰撞 中 高 密 核 物 质 的 状态 方程 .3. 
Burner (CRAB)I 模 型 的 输入 量 构建 关联 函数 : 
500 T T T 
m I J dfx d xS Ga p). Go. py) lp(q, rn 


= - CMF PT1 EoS 

| = CMF PT2 EoS i 
eH CMF PT3 EoS 

= Hard EoS / 
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图 1 (在 线 彩 图 ) 不 同 EoS 给 出 的 密度 依赖 的 平均 场 势 VY (E 

a) 和 相应 的 压强 p (图 b) 
中 四 3435,45]， 通 过 基于 CMF 模型 得 到 的 平均 每 重子 能 
量 计算 出 粒子 坐标 和 动量 时 空 滨 化 所 需要 的 平均 场 势 
能 (外 ， 通 过 改变 势能 极 小 值 点 所 处 的 密度 区 可 以 引入 
在 不 同 密度 区 发 生 强 子 相 与 QGP 相 的 一 阶 相 变 ”435]。 
图 1(a) 展示 了 不 同 EoS 给 出 的 平均 场 势 随 重子 密度 
的 变化 情况 。 由 于 根据 压强 来 讨论 核 物质 状态 方程 的 属 
性 更 加 具有 指导 意义 且 清 晰 明了 ， 因 此 图 1(b) 展示 了 不 
同 EoS 下 压强 随 重子 密度 的 变化 ， 其 中 压强 计算 公式 
Jg [35]; 


7 


Db 


/ 
Pg. T) = Pty T) | ag! 


TALA (3) 


这 里 Pp, T) 是 强 子 的 理想 费 米 气体 压强 ，U(p,) 是 密 
度 依赖 的 单 粒 子 势 能 。 图 1 中 ， 标 准 的 CMF Eos 含有 从 
强 子 气 到 QGP 的 平滑 过 渡 相 变 ，CMF_PT EoSs 中 含有 
从 强 子 气 到 QGP 的 一 阶 相 变 ， 其 中 CMF_PT1 Eos 中 的 
QCD 一 阶 相 变 发 生 在 低 重 子 密度 区 ，CMF_PT3 EoS 中 
的 QCD 一 阶 相 变 发 生 在 高 重子 密度 区 ，CMF_PT2 EoS 
中 的 QCD 一 阶 相 变 发 生 的 密度 区 位 于 二 者 之 间 。 在 参 
考 文献 565] 中 ,作者 探究 了 不 同 Eos 对 集体 流 的 影响 , 通 
过 与 实验 数据 的 比较 排除 了 在 4p, 以 下 发 生 强 相 互 作用 
日 变 的 可 能 性 ， 即 排除 了 CMF. PTI EoS。 因 此 ， 接 下 来 
我 们 将 重点 研究 CMF、CMF_PT2、CMF_PT3、Hard 和 
Soft EoSs 对 两 x 介子 HBT 关联 的 影响 。 


2.2 两 x 介子 HBT 关联 

本 文 将 基于 UrQMD 模型 ， 通 过 以 下 三 步 利 用 两 x 
介子 HBT 关联 抽取 表征 粒子 发 射 源 时 空 属 性 的 HBT 半 
径 及 参数 。 首先, 利用 UrQMD 模型 获得 冻结 点 的 x 介子 


> 


C(k,q)=1+ 
| d^x, S(x;, p) d^x4S (x4, P2) 


这 里 q =p, 一 py Hk = (p, + p,)/2 分 别 是 两 粒子 相对 动 
量 和 平均 动量 。S(x,p) 表示 从 时 空 点 x = (r, A 处 发 射 动 
量 为 p 的 粒子 的 概率 ， 它 包含 了 发 射 源 的 时 空 信息 。 两 
粒子 的 相对 波 函 数 (qr). 包含 粒子 间 的 末 态 相互 作用 
信息 以 及 量子 统计 信息 。 

最 后 ， 采 用 三 维 高 斯 函数 拟 合 构建 的 关联 函数 ， 提 

取出 粒子 发 射 源 的 三 维 半径 : 
C(gqi, qo» ds) =N[1 + Aexp(—R? q? — Roa 
— Ria; — 2R01.49091)]， 
HE N 为 归 一 化 因子 ，4 是 相干 因子 。 这 里 我 们 采用 
纵向 协 变 运动 坐标 系 (the longitudinal comoving system, 
LCMS) 进行 三 维 分 析 。 粒 子 对 的 相对 动量 4 被 分 解 为 
Out, Side, Long 三 个 方向 , 其 中 纵 轴 Long O) 沿 着 束 流 方 
向 , Out(O) 轴 与 粒子 对 的 横向 动量 方向 kr = (Pri + p72)/2 
平行 ，Side(S) #15 Long 轴 和 Out MEH. 

当 研究 对 象 为 带电 粒子 时 ， 构 建 关 联 函 数 的 过 程 中 
必须 考虑 末 态 相互 作用 中 的 库仑 相互 作用 ， 同 时 在 拟 合 
过 程 中 也 应 当 进 一 步 考虑 库仑 修正 347481, AMER 
究 目标 为 x 介子 , 因此 为 了 和 实验 数据 进行 严格 的 对 
进而 抽取 更 为 可 靠 的 核 物 质 属性 信息 ， 在 分 析 过 程 中 我 
们 采用 和 实验 分 析 包 -291 一 样 的 Bowler-Sinyukov 方法 来 
考虑 库仑 相互 作用 的 修正 外， 则 式 (5) 变 为 : 


C(q,, qo: qs) =NI[(l 一 À) + AKcGiny> Ry JA + exp( - R2 q? 
— R? q — Rig? - 2R2 494,))] 


G) 


SES 


sy 


(6) 
不 变相 对 


其 中 Ke 为 qu, 和 Riny 依赖 的 库仑 修正 因子 ， 
动量 quu, 定义 为 

du, = 5 V i - P2}? - (E, ~ E)? 0 
通过 拟 合 得 到 的 R, Ro 和 Rs 分 别 对 应 三 个 方向 的 HBT 
半径 ; Ro 为 交 久 项， 通常 在 分 析 对 称 碰撞 体系 中 大 快 
度 区 间 的 粒子 关联 时 发 挥 作用 。 
3 结果 与 分 析 
3.1 不 含 相 变 的 Eos 的 影响 


图 2 展示 了 碰撞 能 量 为 VSNN = 4.5 GeV 的 Au+Au 
中 心 (0-10%) 碰撞 中 z-z- 的 三 维 HBT 关联 函数 ， 理 论 


相 空 间 信息 ， 这 里 的 冻结 点 指 的 是 动力 学 冻结 点 ， 即 最 
后 一 次 发 生 相 互 作 用 (碰撞 或 者 衰变 ) 的 时 空 点 。 然 后 ， 
将 粒子 在 冻结 点 的 坐标 和 四 动量 作为 Correlation After- 


模拟 计算 均 采 用 了 Hard EoS 模式 ， 红 色 圆 点 为 取 自 文 


"https://web.pa.msu.edu/people/pratts/freecodes/crab/home.html 
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到 2 (ERKE) 碰撞 能 和 


STAR 合作 组 外 


相互 作用 的 影响 ， 即 在 


献 29] 的 实验 数据 。 图 (a-c) 主要 展示 了 粒子 冻 出 前 库仑 


为 VSNN = 4.5 GeV 的 Au+Au 中 心 (0-10%) 碰撞 


NN 


> aT 


q, (MeV/c) 


介子 对 的 HBT 关联 函数 。 实 验 数据 来 自 于 


日 的 实验 数据 由 红色 五 角 星 表示 。 首 先 , 相 较 于 Cascade 


CRAB 分 析 中 不 考虑 末 态 相互 作 


用 中 的 库仑 修正 , 而 在 UrQMD 模拟 过 程 中 分 别 考 虑 ( 黑 


色 方块 ) 和 不 考虑 GEZE) 库仑 相互 作用 ， 可 以 看 


出 冻 出 前 库仑 相互 作用 
库仑 相互 作用 ) 对 xx 


ROSAS! 的 结论 一 致 。 图 (d-f) 展示 的 是 在 UrQMD 模拟 


( 即 UrQMD 模型 模拟 过 程 中 的 
- 关联 函数 的 影响 很 弱 ， 这 与 文 


过 程 中 考虑 库仑 相互 作 
( 蓝 色 三 角 ) 和 不 考虑 ( 


修正 的 结果 ， 可 以 看 到 在 考虑 了 未 态 相 互 作用 中 的 
相互 作用 后 ，x-x- 的 HBT 关联 函数 能 够 更 好 地 描述 
验 数据 。 此 外 , 图 中 给 出 了 一 维 拟 合 结果 ( 蓝 色 实 线 ), 拟 


合 得 到 的 Adj. R-Square 


示 在 对 应 的 图 中 ， 可 以 看 出 关联 函数 的 计算 值 和 实验 值 


], MÆ CRAB 分 析 中 分 别 考 
黑色 方块 ) 未 态 相 互 作用 中 


In IY 
T 


(R2) 值 和 三 维 HBT 半径 分 别 展 


均 能 通过 一 维 高 斯 形式 得 到 较为 完美 的 拟 合 。 因 此 ， 上 
述 关 联 函 数 的 构建 和 拟 合 方法 可 以 用 来 提取 表征 粒子 发 
射 产 时 空 属 性 的 HBT 半径 。 


图 3 展示 了 碰撞 能 量 为 VSNN = 2.4 GeV 下 Au+Anu 


中 心 (0-10%) 碰撞 中 横 动 量 ky 依赖 的 HBT 半径 。 为 了 


阐明 势 修正 及 EoS 对 x 
算 过 程 中 分 别 采用 了 不 


介子 发 射 源 时 空 属 性 的 影响 ， 计 
含 势 修正 的 Cascade 模式 和 含有 


LH 
H 


粒子 间 的 时 间 - 空 间 关 联 ， 导 致 大 kr 区 间 内 的 Ro 被 压 


BS 


模式 ， 考 虑 势 修正 后 的 平均 场 模 式 下 的 结果 能 够 较 好 地 


述 实验 数据 。 其 次 ， 由 于 平均 场 势 修 正 提供 了 更 强 的 
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MA OM BRAN BD o 
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] (al) 
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EX HEKK, £ 


色 点 线 分 别 表示 ，HADES 实验 
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3 (在 线 彩 图 ) VSNN = 2.4 GeV 下 Au+Au Ò (0-10%) fif 


横 动 量 


ky 依赖 的 HBT 半径 Ro, Rs 和 Ruo KIRZA 


F HADES 合作 组 CI 
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24 EREE) 势 修正 对 通过 x- 介子 发 射 源 的 HBT 半 5 “(在 线 彩 图 ) Au--Au 中 心 碰撞 中 x- 介子 发 射 源 的 RO/Rs 
Ro. Ry 和 Rs 的 影响 及 其 随 碰 接 能 量 的 变化 之 比 随 碰撞 能 量 的 变化 。 其 中 实验 数据 取 自 参考 文 
: [28-29] 

K, ^, 区 间 的 Rs 被 抬升 。 最 后 ， 随 着 Eos 硬度 的 d 
增加 (不 可 压缩 系数 Ko 数值 的 增 大 ) ， 密 度 依赖 的 强 。 ”并 且 能 够 较 好 地 描述 现 有 的 实验 数据 ， 同 时 二 者 之 间 的 
子 相互 作 用 带 来 更 大 的 压强 并 进一步 增强 粒子 之 间 的 时 ”差异 随 着 碰撞 能 量 的 升 高 逐渐 减 小 ， 这 意味 着 x 介子 的 
间 - 空 间 的 关联 性 ， 大 ON) kr 区 间 内 的 RoCRS) 被 进 一 HBT 半径 参数 敏感 于 核 物 质 状态 方程 ， 但 随 着 碰撞 能 量 
步 压低 (抬升 )。 基 于 上 述 分 析 ， 可 以 看 出 x 介子 发 射 源 ”的 升 高 ， 强 子 势 的 效应 逐渐 减弱 。 


Ds 


= 


a 


= 的 HBT 半径 对 Eos 比较 敏感 ， 因 此 可 以 用 来 约束 高 密 其 次 ， 对比 Hard EoS 模式 和 CMF EoSs 模式 下 的 结 
区 的 EoS。 R, 发 现在 该 能 区 Hard EoS, CMF EoS 和 CMF_PT3 EoS 
此 外 ， 通 过 束 流 能 量 扫 描 ， 我 们 还 发 现 密度 依赖 的 ”给 出 了 非常 相似 的 结果 。 这 是 因为 在 所 研究 的 碰撞 能 量 


CZ 。 势 修正 的 效应 随 着 束 流 能 量 的 升 高 而 逐渐 降低 。 为 了 清 EEA, 0-10% 的 Au--Au 中 心 碰撞 中 的 重子 密度 范围 大 
sn — 晰 地 展示 势 修正 对 x 介子 发 射 源 HBT 半径 的 影响 及 其 约 为 2 5p 结合 图 1， 可 以 发 现 当 m < 3po 时 ，CMF、 
随 束 流 能 量 的 变化 ， 图 4 展示 了 不 同 能 量 下 Cascade # ^ CMF PT2 和 CMF PT3 EoSs 给 出 了 和 Hard EoS 几乎 相 
RÄ Hard EoS 模式 的 HBT 半径 之 比 。 左 图 为 大 横 动量 。 同 的 压强 ， 因 此 在 VSNN = 2.4 GeV 时 ,除了 Cascade 模 
(275 < kp < 375 MeV/c) 区 间 Ro 的 比值 随 碰撞 能 量 的 。” 式 外 ， 其 余 含 有 势 修正 的 模式 给 出 了 几乎 一 样 的 Ro/Rs; 
变化 ， 中 间 图 和 右 图 分 别 为 小 横 动 量 (75 < kp « 175 ”而 在 重子 密度 小 于 5p KIRAN, CMF 和 CMF_PT3 EoSs 
MeV/c) 区 间 R,. R, 的 比值 随 磁 撞 能 量 的 变化 。 可 以 ”给 出 了 和 Hard Eos 几乎 相同 的 压强 , 因此 在 所 研究 的 能 
看 到 Ro 的 比值 大 于 1 并 且 随 着 束 流 能 量 的 升 高 而 下 降 ，” 量 范围 内 三 者 给 出 了 相似 的 Ro/Rs 随 碰撞 能 量 的 变化 。 
R, 和 Rs 的 比值 随 着 碰撞 能 量 的 升 高 而 上 升 ， 即 考虑 势 最 后 ， 对 比 不 同 CMF EoSs 的 结果 ， 发 现 随 着 碰撞 
修正 后 的 大 kr 区 间 内 的 Ro 和 小 ky 区 间 的 Rp. Ry 在 。 ”能量 的 增 大 ， 相 比 于 CMF 和 CMF_PT3 EoSs 模式 的 结 
数值 上 分 别 小 于 或 大 于 Cascade 模式 的 结果 ; 同时 随 着 R, HF CMF_PT2 EoS 中 强 子 气 到 QGP 相 的 一 阶 相 变 
束 流 能 量 的 升 高 ， 势 修正 的 效应 逐渐 减弱 ， 两 种 模式 下 ” 发 生 在 相对 较 低 的 密度 区 间 , CMF_PT2 EoS 模式 下 的 结 
结果 之 间 的 差异 也 逐渐 减 小 ， 这 与 参考 文献 59 中 正 反 果 随 碰撞 能 量 的 升 高 逐 潮 增 大 并 最 终 靠近 Cascade 模式 
粒子 椭圆 流 之 间 的 流 差 随 着 束 流 能 量 的 升 高 而 逐渐 降低 。 ”的 结果 。 这 里 需要 特别 指出 的 是 ，Cascade 模式 可 以 看 做 
相 一 致 。 一 种 非常 软 的 Eos 模式 ， 即 发 生 了 相 变 的 模式 。 因 此 在 
相对 较 高 的 能 量 下 Cascade 的 结果 与 CMF_PT2 的 结果 

3.2 仿 相 变 的 Eos 相似 

为 了 更 加 清晰 的 展示 势 修 正 对 两 x 介子 HBT 关联 综 上 所 述 ， 可 以 得 出 结论 : 直到 4-5po XT, a jT 
及 工 介 子 发 射 源 时 空 属 性 的 影响 及 其 随 碰 撞 能 量 的 变 — 子 发 射 源 的 HBT 半径 参数 Ru/Rs 对 EoS 非常 敏感 ; 现 
化 , 图 5 比较 了 不 同 EoS F, a7 介子 发 射 源 的 半径 Ro/Rs ”有 的 HBT 实验 数据 支持 一 个 在 低 密度 时 表现 为 相对 较 
之 比 随 磁 撞 能 量 的 变化 ， 并 且 与 实验 数据 进行 了 对 比 硬 然 后 在 高 密 区 由 于 相 变 而 逐渐 软化 的 EoS; 同时 在 所 
首先 ， 通 过 对 比 Cascade 模式 和 Hard EoS 模式 下 ”研究 的 能 量 范围 内 ，CMF_PT3 EoS 中 强 子 气 到 QGP fH 
的 结果 ， 发 现 引 入 强 子 相互 作用 势 后 ，Ro/Rs 的 值 降 低 ， ”之 间 一 阶 相 变 发 生 所 需要 的 高 密 条 件 从 未 真正 达到 |。 
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图 6 (在 线 彩 图 ) 不 同 状态 方程 下 ，Au+Au 中 心 碰撞 (0-109) ”图 7 (在 线 彩 图 ) Au+Au 中 心 碰撞 (b=0 fm) 中 ， 图 (a): 通过 
中 a7 介子 的 平均 冻 出 时 间 (0). 随 碰 撞 能 量 的 变化 弦 激 发 和 碎 裂 、 共 振 衰变 产生 的 x 介子 数目 占 总 x 介子 
、 数目 的 百分比 ; 图 (b): 重子 -介子 多 重 数 比 值 随 碰撞 能 
3.3 工 介 子 冻 出 时 间 分 布 最 
量 的 变化 
为 了 理解 不 同 Eos 对 两 x 介子 HBT 关联 及 x 介子 。 时 间 将 会 增加 B11。 这 里 需要 指出 ， 在 高 能 量 重 离子 磁 


发 射 源 HBT 半径 的 影响 ， 接 下 来 将 对 不 同 EoS 下 x 介 T, x 介子 的 产生 及 演化 过 程 非常 复杂 , 一 方面 纺 激 发 和 
子 的 冻 出 时 间 进 行 系统 地 分 析 。 图 6 给 出 了 Au+Au Pò 。 碎 裂 可 以 直接 产生 x 介子 并 被 直接 冻 出 或 者 通过 M M. 
碰撞 (0-10%) 中 x- 介子 的 平均 冻 出 时 间 随 碰撞 能 量 的 。 MB 相互 作用 重新 产生 弦 并 进一步 碎 裂 ， 同 时 弦 激 发 和 
变化 。 首 先 ， 可 以 看 到 在 低 碰 撞 能 量 (VSNN = 24 GeV) 碎 裂 也 可 以 产生 共振 态 并 通过 共振 衰变 产生 x 介子 或 者 
时 ， 考 虑 势 修正 后 (Hard, CMF, CMF_PT2 和 CMF PT3) ”通过 BB 相互 作用 重新 产生 弦 并 进一步 碎 错 。 男 一 方面 
的 x- 介子 的 平均 冻 出 时 间 基 本 一 致 ， EF Cascade 共振 衰变 产生 的 x 介子 和 共振 态 也 可 能 经 历 上 述 过 程 。 
模式 下 x- 介子 的 平均 冻 出 时 间 ， 这 是 因为 在 该 能 量 下 ”第 二 ， 自 由 度 的 变化 。 如 图 7(b) 所 示 ， 在 考虑 和 不 考虑 
CMF, CMF_PT2 和 CMF_PT3 EoSs 给 出 了 和 Hard EoS 。” 势 相互 作用 的 两 种 模式 下 , Au+Au FEUCMIÉTRE (b=0 fm) 中 


相似 的 密度 依赖 的 排斥 的 势 相 互 作用 ， 而 硬 的 密度 依赖 。 重子 -介子 多 重 数 比值 随 碰撞 能 量 的 变化 表现 一 致 ,并 

的 势 相 互 作用 将 提供 较 强 的 压强 (排斥) 并 带 来 更 强 的 。 在 VSNN RA GeV 时 Nparyon X Nmesono 因此 , 在 介子 

膨胀 ,因此 x- 介子 被 更 早 的 冻 出 。 其 次 ， 随 着 碰撞 能 量 。 ”的 平均 冻 出 时 间 和 粒子 发 射 源 体积 的 激发 函数 上 所 体现 
-— 的 升 高 ,在 CMF_PT2 EoS 模式 下 , 由 于 相 变 的 发 生 , 状态 。 ”出 来 的 最 小 值 可 能 由 不 同 的 粒子 产生 机 制 之 间 的 苋 争 与 
号。 方程 逐渐 软化 ， 因 此 其 结果 逐渐 靠近 Cascade 模式 下 的 。 ”改变 或 重子 -重子 相互 作用 与 介子 (重子 ) -介子 相互 作 


结果 。 此 外 , 在 VSNN =7.7GeV I, CMF 和 CMF_PT3 ”用 的 竞争 与 转变 来 解释 。 
EoSs 的 结果 位 于 Hard EoS 和 CMF. PT2 EoS 的 结果 之 

间 ， 这 是 因为 在 该 碰撞 能 量 下 CMF 和 CMF_PT3EoSs ^ 总 结 与 展望 

比 Hard EoS 软 ， 但 比 CMF_PT2 EoS 硬 ， 如 图 1 中 所 示 。 基于 极端 相对 论 量子 分 子 动力 学 模型 ， 本 文 讨论 了 
有 意思 的 是 ， 在 图 6 中 无 论 是 否 考虑 强 子 势 或 者 相 ”密度 依赖 的 核 物质 状态 方程 和 相 变 对 两 x 介子 HBT X 
变 , xz” 介子 的 平均 冻 出 时 间 (0) 随 能 量 的 变化 在 VSNN mo 联 和 x 介子 发 射 源 的 HBT 半径 参数 的 影响 。 主 要 结论 
4 GeV 时 均 出 现 了 极 小 值 点 。 同 时 ， 我 们 注意 到 在 现 如 下 : 首先 ， 厘清 了 zz 介子 冻 出 前 后 库仑 相互 作用 对 x 


有 的 HBT 数据 分 析 中 ，z 介子 发 射 源 的 体积 VV, = MTHAI HBT 关联 函数 以 及 HBT 半径 的 影响 ， 发 现 冻 
(22)? R? Rion 的 激发 函数 也 在 VSNN m 4 GeV 呈现 出 前 的 库仑 相互 作用 对 sa 的 HBT 半径 和 Ro/Rs 之 
出 了 相似 结构 ?3 我们 提出 以 下 两 种 原因 来 解释 这 一 ” 比 有 着 很 弱 的 影响 ， 而 末 态 库仑 相互 作用 却 对 x-x- 的 


极 小 值 迷 团 。 第 一 ， 粒 子 产 生机 制 的 转变 。 如 图 7(a) 所 ”HBT 关联 函数 有 着 较 强 的 影响 。 其 次 , 通过 对 比 Cascade, 
示 ， 在 UrQMD 模型 中 ， 当 束 流 能 量 [su S 4 GeV IN, Hard EoS 和 Soft EoS 模式 下 的 结果 发 现 密度 依赖 的 排斥 
z 介子 的 产生 机 制 主要 以 共振 衰变 为 主 ， 随 着 能 量 的 升 ”的 势 修正 会 进一步 增强 粒子 之 间 的 时 空 关 联 性 ， 导 致 大 
高 ， 通 过 弦 激 发 和 碎 裂 产生 的 a 介子 逐渐 增多 ， 并 且 更 ok, 区 间 的 Ro 降低 ， 小 kr 区 间 的 Rs 上 升 ， 从 而 导致 
多 的 x 介子 将 在 经 历 过 再 散射 后 被 冻 出 ， 因 此 平均 冻 出 。” Ro/Rg 之 比 降低 并 且 能 够 较 好 地 描述 实验 数据 。 最 后 ， 
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在 4~5 倍 饱和 密度 以 下 ， HBT 半径 和 Ro/R, 之 比 敏感 


硬 ， 


f, 


旧 存 在 很 大 的 不 确定 性 。 因 此 ， 
实验 数据 , 相 


于 核 物质 状态 方程 ， 通 过 与 实验 数据 的 比较 ， 发 现 现 有 
的 HBT 实验 数据 支持 一 个 在 相 


7 


x 


对 低 密 区 (< 4p9) 表现 偏 
然后 在 高 密 区 由 于 相 变 而 逐渐 软化 的 核 物质 状态 方 
并 且 排 除了 在 < 4p 以 下 发 生 一 阶 相 变 的 可 能 性 。 

但 是 , 在 本 文 所 研究 的 磁 撞 能 量 范 围 内 , 实验 数据 依 
方面 需要 更 加 精确 的 
信 未 来 德国 GSI 上 FAIR 的 CBM 实验 和 俄 


罗斯 杜 布 纳 NICA 上 的 MPD 实验 以 及 美国 RHIC-STAR 


固定 革 项 目 将 提供 


i 加 精确 的 实验 数据 ; 另 一 方面 也 需 


HE 


要 更 多 的 理论 研究 来 探讨 模型 等 不 确定 性 ， 如 初始 化 结 


构 效 应 、 量 子 效应 、 粒 子 势 相 
子 关联 效应 所 带 来 的 影响 ， 其 中 介子 (K, r-T AH 


互 作用 及 动力 学 效应 对 粒 


互 作 


对 两 粒子 HBT 关联 的 影响 的 研究 正在 进行 中 。 
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CNPC2023:The High-density Equation of State in Heavy-ion Collisions: 
Constraints from Two-pion HBT Correlation 


LI Pengcheng'?, WANG Yongjia!, LI Qingfeng!, ZHANG Hongfei? 


(4. School of Science, Huzhou University, Huzhou 313000, Zhejiang, China; 
2. School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, Gansu, China; 
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Abstract: The investigation of the equation of state (EoS) of nuclear matter, especially at high baryon densities is 
one of the hot topics in the frontier of nuclear physics. The impact of the EoS at 2 ~ 5 times saturation density po 
on the two-particle correlation is discussed with the ultra-relativistic quantum molecular dynamics (UrQMD) model. 
Focusing on the two z Hanbury-Brown-Twiss (HBT) correlations, by adopting different EoSs, the effects of potential 
interaction and phase transition on the HBT correlation and the spatiotemporal properties of the emission source of z 
are investigated. The results show that below 4 ~ 5p,, the HBT radius and parameters are sensitive to the stiffness of 
the EoS. By comparing with the experiment data, a phase transition with a significant softening of the equation of state 
below 4 times nuclear saturation density can be excluded using HBT data, and the available data on the HBT radii in the 
investigated energy region favor a relatively stiff EoS at low densities, which then turns into a soft EoS at high densities. 
These results highlight that the pion's HBT radius and parameters are sensitive to the stiffness of the equation of state, 
and can be used to constrain and understand the equation of state in the high baryon density region. 


Key words: heavy-ion collisions; equation of state; transport model; HBT interferometry 
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